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ABSTRACT 

We derive the conditions for a backflow toward the central star(s) of circumstellar material to 
occur during the post-asymptotic giant branch (AGB) phase. The backflowing material may be 
accreted by the post- AGB star and/or its companion, if such exists. Such a backflow may play a 
significant role in shaping the descendant planetary nebula, by, among other things, slowing down 
the post- AGB evolution, and by forming an accretion disk which may blow two jets. We consider 
three forces acting on a slowly moving mass element: the gravity of the central system, radiation 
pressure, and fast wind ram pressure. We find that for a significant backflow to occur, a slow dense 
flow should exsist, such that the relation between the total mass in the slow flow, Mj, and the solid 
angle it covers Q, is given by {Mi/f5) ^ O.IMq, where /3 = r2/47r. The requirement for both high 
mass loss rate per unit solid angle and a very slow wind, such that it can be decelerated and flow 
back, probably requires close binary interaction. 

Key words: Planetary nebulae: general — stars: AGB and post-AGB — stars: mass loss — circum- 
stellar matter 
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1. INTRODUCTION 

More and more supporting observations (Sahai k, Trauger 1998; Kwok, Su, & Hrivnak 1998; 
Hrivnak, Kwok, & Su 1999; Kwok, Hrivnak, & Su 2000; Huggins et al. 2001) and theoretical 
considerations (Soker 1990; Soker 2001) are accumulated in support of the view that significant 
shaping of the circumstellar material takes place just before, after, and during the transition from 
the asymptotic giant branch (AGB) to the planetary nebula (PN) phase. Both the wind and 
radiation properties are significantly changed during these stages. As the star is about to leave 
the AGB the mass loss rate increases substantially, up to ~ 10~^ — W'^Mq yv"^. This wind was 
termed superwind by Renzini (1981); because of confusion with other winds termed superwinds, 
hereafter we'll refer to this wind as FIW, for final intensive wind. After the star leaves the AGB 
the mass loss rate decreases down to 10^^ Mr.) yr~^, and its velocity increases from ^ 10 km s^^ 
to fewxlO^ km s~^ at the PN phase. Simultaneously, the effective temperature increases, and the 
post-AGB star starts to ionize the nebula around it, when by definition the PN phase starts. The 
changes in the mass loss rate and velocity of the wind are accompanied by a change in the wind 
geometry, e.g., jets and bipolar structures are formed (Soker 1990; Sahai Sz Trauger 1998; Kwok 
et al. 2000). The change in the wind geometry may result from an intrinsic processes in the AGB 
and post-AGB mass-losing star, or from processes in an accreting companion (e.g., Soker 2001). 
The increase in the wind velocity leads to collision of winds (Kwok, Purton, &; Fitzgerald 1978), 
which can lead to instabilities during the post-AGB (proto-PN) stage (Dwarkadas Sz Balick 1998). 
As the central star starts ionizing the nebula, an ionization front propagates outward, and plays 
a significant role in shaping the nebula, both in the radial direction (e.g., Mellema &; Prank 1995; 
Chevalier 1997; Schonberner &; Steffen 2000), and in the transverse directions (e.g., Mellema 1995; 
Soker 2000b). 

In the present paper we examine yet another process which may occur during the transition 
from the AGB to the PN stage, i.e., after the FIW ceases and before ionization starts. This is a 
backflow of a fraction of the dense wind toward the central star(s). This backflowing material may 
be accreted by the central star and/or its companion. The idea of accreting backflowing material 
during the post-AGB phase was raised before to explain and account for: a possible mechanism for 
the formation of jets from an accretion disk (Bujarrabal, Alcolea, & Neri 1998); a slower evolution 
along the post-AGB track (Zijlstra et a;. 2000); and post-AGB stars depleted of refractory elements 
which compose the dust particles (e.g.. Van Winckel et al. 1998) by accretion of a dust-depleted 
circumstellar gas (Waters, Trams & Waelkens 1992), most likely in binary systems (Van Winckel 
1999). 

The goal of the present paper is to explore the conditions required for a backflow to occur such 
that it plays a non-negligible role in the post-AGB evolution. The conditions are derived in §2, 
while the implications for the processes mentioned above, as well as other processes, together with 
a short summary, are in §3. 
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2. CONDITIONS FOR A BACKFLOW 



In this section we derive the conditions for a backfiow to occur during the post-AGB phase. 
We assume that a very slow flow exists along some directions, e.g., in the equatorial plane, and 
consider the conditions for some of this material to flow inward and be accreted by the central 
system, before ionization starts. After being ionized, any dense cool gas will expand and will be 
pushed outward by radiation and ram pressure (see below). We do not consider the deceleration 
of the slowly outward moving mass element, but simply assume that if there is a slowly outward 
moving mass element, it will reach zero radial velocity at some radius. 

We therefore consider a mass-element Mj with a constant density pi within a solid angle f2 
and a radial extention Ar at a distance r :§> Ar from the central star, such that 

Mi = piQr'^Ar. (1) 

The mass element can also be a shell where Q = 47r. If the sound crossing time Ar/cg, where Cg 
is the sound speed, is shorter than any other time scale in the process, we can take the shell to 
move more or less coherently. Three relevant forces are acting on the mass element in the radial 
direction. The gravitational force of the central star(s) 

where M is the total mass of the central system, a binary system or a single star. The fast wind 
blown by the central star and its radiation push outward. The force due to the fast wind's ram 
pressure, assuming it is much faster than the slow wind velocity, is given by 

fw = Pw{r)vlnr'^ = MwVwP, (3) 

where Pw{r) = M«,/(47rr^t;tj,), and Mu, are the density, velocity and mass loss rate (defined 
positively) of the fast wind, and (3 = ^l/An. The radiation imparts a force of 

f, = ^P{l-e-^), (4) 

c 

where is the luminosity of the central system, c the speed of light, and 

^ = ^^"^^ = ^ (oollo) do cm"^ g-i ) (ioTau) 

is the optical depth of the mass element, and we used equation (1) for the density. We scale the 
opacity K{r) with a typical value for AGB stars (Jura 1986; Winters et al. 2000), and assume 
that the fast wind inner to the mass element absorbs a negligible fraction of the radiation. For 
convenience we define three dimensionless variables. The ratio of maximum radiation pressure to 
the fast wind ram pressure 

^ My,vy, UOOOL©; i IO-6M0 yr-1 i boO km s" V ' ^' 
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which does not depend on r. We scale the fast wind mass loss and velocity as appropriate for a 
post-AGB star before it starts ionizing the nebula. The ratio of the force due to radiation and wind 
to that of gravity depends on r, both through the dependence of optical depth x cind gravity on r. 
We define t] to be this ratio at a scaling radius ro 

where 

fc(ro) =/?[! + (/(I -e-^)]. (8) 
The equation of motion for the mass element can be written as 

(9) 




We take ro to be the radius at which the radial velocity of the mass element is zero. When there 
is only gravity, the free-fall time from r = ro to the center r = 0, with v{ro) = 0, is 

We define the dimensionless variables 

T = t/tff and x = r/ro, (11) 



and write the equation of motion (9) in the form 

= ^{ri-x-^). (12) 



Assuming that rj is constant and does not depend on x allows us to integrate once the last equation 
to give the velocity 

^^^ = -lirix + x-'-l-r))y\ (13) 

where we substituted the initial condition t;(l) = 0. This can be integrated analytically for constant 
values of r/ = 1 and rj = 0. The case rj = gives the free fall solution. The time left for the object 
to fall from x to a; = is given by 

Tf{x) = ^ (sin-1 a;i/2 - [x{l - x)]^/^) . (14) 

As expected from our scaling t(1) = 1, i.e., the free fall time from r = tq to r = is tff. For r] = 1 
and with ^;(1) = 0, the fall time from x = 1 is infinite. This is because the wind and radiation 
outward-acceleration equals the gravity inward-acceleration. However, the time left to fall from x 
to a; = 0, with v{l) = 0, is finite, and is given by 

r,W = ^(ln[(l + x'/=)/(l -.■/.)-.,_ 2,./.). (15) 
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Comparing the last two equations, we find as expected, that radiation pressure and wind's ram 
pressure slow down the inflow. For example, for i] = 0, i.e., no wind and radiation pressure, the 
time left to fall from x = 0.9 to .t = is r = 0.604, while for r/ = 1 it is r = 1.11. In both cases 
the initial condition is v{l) = 0. These cases are less interesting, since for i] = 1 the flow actually 
stagnates at x = 1. More interesting is the case of v{l) = and < rj < 1, since for 77 > 1 it will 
be push away from the center. We numerically integrated equation (13) for these conditions and 
for constant values of rj. We find the backfiow time to be t/(1) = 1.28, 1.52, and 2.08, for 77 = 0.5, 
0.7, and 0.9, respectively. For r] = 0.8, for example, the fall back times from Xi = 1, 0.9 and 0.8 
are Tf = 1.72, 1.23 and 0.925, respectively, where v{xi) = in these cases. For ij = the backflow 
times for the same initial conditions are 1, 0.85, and 0.72, respectively. The conclusion from the 
numerical values cited above is that the typical backflow time from r ~ Tq is the free fall time ~ t// 
at ro, but because of the radiation and wind pressures the region from which this is the backflow 
time is much larger than the r/ = cases, extending from ro down to ~ 0.7 — 0.8ro for 77 ;^ 0.8. So 
the question is what is the value of rg for which r] = 1. Below and close to this radius the gas falls 
back in a time ^ tff, while it is accelerated away for larger radii. From equation (7) we find 

M \ / Mi \ / \ / Vy, \-V k ^ 



ro = —. = 430 AU 

V kMwVw J 



IMq/ VO.OIM0; 110-6Mq y VlOOkms-V VO.I 



(16) 



For the backflowing mass to influence the evolution signiflcantly, we required the fall back time 
to be ^ 1000 yr. For the typical parameters used in equations (5) and (6) we flnd from equation 
(8) that k{r) ~ 2/3; using the time scale given by equation (10) in equation (16) gives the desired 
condition 

-J O.lMc (17) 

where as before, /3 = Q/Att, and Q is the solid angle covered by the dense backflowing material. 
The last condition is limited by a maximum density, since a large value of the backflowing mass 
Mj and small value for /? means a very high density. We now estimate a reasonable value for the 
density. We assume a very slow equatorial flow, with a speed of ~ 1 km and with a mass loss 
rate per unit solid angle of rhs = Mg/^T^- The density of the mass elements formed by this wind is 

The minimum density is that for which the fast wind compresses the dense cool wind such that the 
ram pressure pwvl^ equals the thermal pressure of the cool gas. For a molecular gas, we find this 
density to be 

«p ^ .4 ^ 10- (^) (^^^) {-^) {^y\ c.-» (19) 

where Cj and Tj are the isothermal sound speed and temperature, respectively, of the cool gas. We 
find that a density of pi ~ 10"^^ g cm~^ is reasonable. Using equation (1) and the definition of /? 



-6- 



gives 

f = iw^w) iwA^y (^) 

Condition (17) is met for the density given by equation (18), but this requires a very high mass loss 
rate per unit soHd angle. If /3 = 0.1 this requires a total mass loss rate of M = 10~'^ Mq yr^^, but 
concentrated in particular directions, probably in the equatorial plane. All these considerations 
strongly suggest an equatorial dense and slow flow, such as expected in a close binary system 
(Mastrodemos Sz Morris 1999; Soker 2000a). A very fast rotation can also form such a wind 
(Bjorkman & Cassinelli 1993), but then a binary companion is needed to substantially spin-up the 
envelope. 



3. IMPLICATIONS AND SUMMARY 



Despite the assumptions and simplifications in deriving condition (17), wc feel that the results 
obtained in the previous section and the implications discussed in this section are quite general. We 
find that for a backflow to occur on a time scale of iacc ^ 10^ yr after the termination of the AGB, 
so that it has a non-negligible role in the post-AGB evolution, the following conditions should be 
met: (1) The total backflowing mass should be larger than the combined mass lost in the wind and 
that burned in the core. For a post-AGB mass loss rate of ~ 1O~^M0 yr~^ the nuclear burning is 
negligible, and the total required mass is Mace — 10^'^(tacc/l5 000 yr)MQ. (2) The backflowing mass 
should have a very low, ~ 1 km s~^, terminal velocity, so that eventually it will be decelerated to 
zero velocity, and flow back. (3) Condition (17) on the ratio of the mass of the mass-element and 
the solid angle it covers /3 = CI /An should be met. (4) For reasonable densities (eq. 18), and the 
required mass, we find (eq. 20) that the mass range is Mj ~ 0.1 — 10~^Mq, and the appropriate 
solid angle covered by the backflowing mass is 1 < ^ ^ 10~^. We can take the typical values to be 
Mi ~ O.OlMo and /3 = 0.1. 

These mass loss properties required that (i) the flow be concentrated along particular direc- 
tions, and (m) have an inefficient acceleration by the stellar radiation. In a previous paper (Soker 
2000a) two mechanisms which lead to such a flow were discussed. In the flrst mechanism proposed 
by Soker (2000a) magnetic cool spots (as in the Sun) are formed on the surface of slowly rotating 
AGB stars. The lower temperature above the spots enhances dust formation, which during the 
final intensive wind (FIW) may lead to an optically thick wind. This means an inefficient radiative 
acceleration, hence a slow flow above the spot. If the spot is small, material from the surrounding 
flows into the shaded region and accelerates the slow flow. However, if the spot is large, material 
from the surroundings will not accelerate the flow much, and it will stay slow. In the present 
paper we find that the condition for the slowly moving material to fiow back is that the spot has 
P ^ 0.01, which for a circular spot means a radius of -Rgpot ^ 0.2i?*, where i?* is the stellar radius. 
The required mass in the slow flow is given by equation (17). 
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In the second mechanism a high density in the equatorial plane is formed by a binary inter- 
action, where the secondary star is close to, but outside the AGB envelope. This mechanism is 
supported by observations, e.g., slowly moving equatorial gas is found around several binary post- 
AGB stars (Van Winckel 1999, and references therein). In the process proposed by Soker (2000a) 
the strong interaction between the two stars forms a dense equatorial outflow which is optically 
thick, leading to an inefficient radiative acceleration and a very slow equatorial flow. For a very 
massive and significant backflow, with a mass of Mj ^ O.OIMq and (3 ^ 0.1, to occur, a binary 
mechanism is required. For the cool spots model to form such a massive slow flow, several large 
spots are required. This means a strong magnetic activity, which probably requires the AGB star 
to be spun-up by a stellar companion. We therefore argue that in both mechanisms a binary com- 
panion is required to cause a massive flow, such that it may last for tacc ^ 10^ yi's, possibly by as 
long as ~ 10^ yrs in extreme cases. 

Such a backflow may have the following effects. If it has a large specific angular momentum, 
as expected in strongly interacting binary systems, the backflowing material may form an accretion 
disk around one or two of the two stars. The disk(s) may blow jets or collimated fast winds (CFWs), 
which will play a significant role in shaping the circumstellar material. Such a possibility was briefly 
mentioned by Bujarrabal et al. (1998) for the proto-PN Ml-92. The accreted mass may slow down 
the post-AGB evolution, as suggested by Zijlstra et al. (2000) for some OH/IR stars. Zijlstra et al. 
(2000) termed these stars retarded stars, and argue for a delay as long as 10^ yrs by accretion from 
a near-stationary reservoir. They bring supporting observations, but don't consider the formation 
of such a reservoir of mass. The results of the present paper put the idea of Zijlstra et al. (2000) on 
a more solid ground. Another effect attributed to backflow accretion is the formation of post-AGB 
stars depleted of refractory elements which compose the dust particles (e.g.. Van Winckel et al. 
1998; Waters, Trams &: Waelkens 1992; Van Winckel 1999). The separation of dust from the gas 
was not considered in the present paper. 

Finally, we speculate on another plausible effect of the dense backflowing gas. The central 
star's wind is shocked when it hits the dense material. If some dense backflowing blobs survive 
long into the PN phase, when the central star's wind velocity is ^ 10^ km s"-*^, then there will be 
a hard X-ray emission from the post-shock fast wind material. The dense blob will be close to the 
central star, making the hard X-ray emitting region hard to resolve. It is not clear if this compact 
hard X-ray emitting region can explain the recent Chandra observations of a "point source" in the 
centers of the Helix (NGC 7293) and Cat's Eye (NGC 6543) PNe (Guerrero et al. 2000). 
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